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DOUBLE ENERGY RADIOGRAPHY METHOD, AND CALIBRATION 
DEVICE FOR THIS METHOD 

The present invention relates to a double energy 
radiography method, together with an appropriate 
calibration device for this method. 

Double energy radiography consists of exposing an 
subject or a patient to radiation with two different 
energies, for which the attenuation properties of the 
materials constituting the subject or the patient do 
not vary in the same way. Using a model of their 
attenuation functions, one thus obtains information 
about the density and nature of the material crossed 
through. In particular, in the case of bone 
densitometry, the bone mass crossed can be calculated 
by distinguishing it from the contribution of soft 
tissues to attenuation of the radiation. 

The theoretical basis of the method will be 
resumed briefly below. 

The flux 0, after crossing a length 1 of a 
material with linear attenuation coefficient p from an 
initial radiation flux 0 O , is equal to 0=0 O e~ ul . The 



where I 0 and I are signals measured by a same detector 
under the fluxes 0 O and 0. In the case of a complex 
subject composed of a large number of materials (of 
index i) , each one contributes to attenuation according 
to its length Li crossed by the rays. 

However, each material can be expressed, for its 
attenuation property, as a linear combination of two 
base materials, according to the formula: 



attenuation measurement is named m, equal to 
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P = kiji! + k 2 p 2 , 
where ki and k 2 are constant coefficients, and pi and P2 
represent the attenuation of these base materials, and 
the equivalent lengths Ai and A 2 of the base materials 
crossed by the radiation by: 

By means of these equivalent lengths, the base 
materials can represent a subject even if in reality 
its composition is much more complex. Even if the 
materials composing the subject are different from the 
base materials, the breakdown has a meaning. In the 
case of examining living beings, the classic base 
materials are Plexiglas (polymethacrylate) to simulate 
soft tissues and hydroxyapatite to simulate bone 
tissues . 

The system of equations linking the measurements 
to the attenuation pi and to the equivalent lengths Ai 
of the base materials is linear and therefore simple to 
solve if the radiation is monochromatic. But this is 
not so for real situations, and the equivalent lengths 
are then given by more complicated mathematical models, 
such as: 

2 2 

A x = a 0 + a x .mBE + a 2 .mHE + a 3 .mBE .mHE + a 4 .m BE +a 5 .m HE 

2 2 

A 2 = b 0 + bi.mBE + fy.mHE + b3.mBE.mHE + b 4 .m BE +b s .m HE 

which generally provide sufficient precision, and where 
mBE and mHE represent the attenuation measurements at 
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high and low energy, and "a" and M b" are coefficients 
which must be calculated beforehand by calibration. 

This calibration requires a device which is often 
called a "phantom" and which is composed of base 
materials as described above and chosen in order to 
simulate the subject to be measured as closely as 
possible. These materials are distributed in the 
phantom in such a way as to provide regions where the 
lengths of materials, crossed by the radiation, are 
distributed differently. 

A classic phantom is built up in the form of a 
staircase for each of the measuring materials, and the 
stairs are superposed in such a way that their steps 
are perpendicular. By making vertical rays cross the 
steps, one thus obtains all the required combinations 
between the diverse thicknesses of the two materials. 

Another phantom is proposed in the US patent 
5 493 601 and comprises a series of tubes converging 
towards the source of the beam and provided with 
heights divided unequally between the two materials. 
The aim is to provide more exact measurements than 
previously, in particular because of the convergence of 
the tubes towards the source, which makes the length 
the rays pass through in the tubes coincide perfectly 
with the total height of the latter, and also reduces 
the scattered radiation by containment in the tubes. 
The measurements are thus a direct expression of the 
attenuation of the radiation through the phantom and 
make it possible to calculate the coefficients required 
reliably. 
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However, the same cannot be said for the 
measurements produced through the subject to undergo 
radiography, for which the favourable conditions 
described above cannot be assembled and the scattered 
radiation must not be neglected; it is even especially 
important for sources with a conical beam, often more 
so than for the primary radiation, which has followed a 
straight line trajectory from the source to the 
detector and is the only one useful for the 
measurement. The estimation and then correction of this 
scattered radiation requires specific processing. Thus, 
amongst other methods, it is possible to measure it 
separately and then to subtract it from the total 
radiation received by the detector. In order to do 
this, one uses a network of absorbent elements, such as 
lead balls arranged in such a way as to form a grid. 
The balls are glued on Plexiglas and set in a 
homogeneous pattern, regular, for example, for the 
lines and columns so as to allow interpolation between 
the balls, between the radiographed subject and the 
detector's network of pixels. The rays intended to pass 
through these balls are intercepted totally, so that 
only the scattered radiation reaches the pixels of the 
detector located in the path of these rays. Digital 
interpolations then make it possible to estimate, with 
sufficient precision, the distribution of the scattered 
radiation over the whole of the detector's pixels. This 
convenient method, however, has the disadvantage of 
having to submit the subject to a second irradiation, 
which may be difficult for living beings to accept. 
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On the other hand, one can also prevent the 
scattered radiation from reaching the detectors by 
providing them with a strict collimation for 
interception. The methods using this technique 
necessitate the use of source scanning which implies a 
long acquisition time and the risk that the subject may 
move during measurement. 

Digital methods of various types have also been 
proposed to correct the influence of scattered 
radiation, but generally they are only suitable for 
particular subjects. It has to be concluded that the 
patent quoted above does not appear to propose any 
notable progress, since the scattered radiation is only 
separated from the measurements at the calibration 
stage . 

The aim of the present invention is to correct the 
influence of the scattered radiation in another way, by 
allowing it to appear during calibration, but in an 
analogous manner to its behaviour during radiography of 
the subject, which makes it possible to correct it by 
an identical digital method, without fear of 
discordance of quality for this method. 

To resume, the invention concerns above all a 
calibration device for a double energy radiography 
system, comprising a range of blocks with different 
thicknesses of a first material, characterised in that 
the blocks are provided with recesses and in that 
furthermore the device comprises inserts to fill up the 
recesses and including different distributions of 
heights, the heights and the thicknesses being 
considered in an identical direction, between the first 
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material and a second material; together with a 
radiography method with a double energy conical beam, 
comprising a thickness estimation of materials of a 
radiographic subject by a digital combination of 
measurements of attenuation of energies, involving a 
coefficient calibration of the combination, 
characterised in that the calibration is measured by a 
radiography of a calibration device conforming to any 
one of the preceding claims, and that scattered 
radiation affecting the radiography of the calibration 
device is estimated while providing an estimation 
criterion used afterwards to estimate a scattered 
radiation affecting the radiography of the subject. 

The invention will now be described with reference 
to the figures: 

- figure 1 shows a flow chart of the invention; 

and 

- figures 2 and 3 show views from the side and 
from above of a calibration phantom used to satisfy the 
requirements described above. 

Reference is now made to figure 1. 

At stage El, a calibration acquisition is carried 
out by a radiography through the phantom described 
below, which gives rough measurements at high and low 
energy. At stage E2, preliminary processing for 
estimating and subtracting the scattered radiation from 
these measurements is made in order to calculate the 
corrected measurements mHe and mBe of the above 
formulae. Next, stage E3 consists of an estimation of 
the parameters of the model, particularly the 
polynomial coefficients a and b, from measurements and 
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the known thicknesses the radiation crosses in the 
phantom. Furthermore, an acquisition, through the 
measured subject, is carried out at stage E4, then a 
preliminary processing is made at stage E5; it is 
similar to that of stage E2 and its aim is also to 
estimate and then subtract the scattered radiation 
which affects the measurements. The final stage E6 
consists of using the corrected measurements made 
through the subject and applying the model, and above 
all the coefficients a and b calculated at stage E3, in 
order to deduce the equivalent lengths of the base 
materials al and a2 in the subject. 

An example of a phantom is shown in figures 2 
and 3 . 

It comprises several blocks of different, 
thicknesses of one of the base materials, in particular 
Plexiglas. An assembly of separate blocks could be 
used, through which the measurements would be made 
successively, but one can also unite the blocks to form 
a solid with a stepped profile 1 comprising different 
layers, four in this case, and numbered from 2 to 5 
from bottom to top, the layers usually finishing by 
dressed tapered faces 6, so as to reduce the high 
frequency components of the scattered radiation. 
Furthermore, rows of inserts are provided under the 
different blocks, filling the recesses made here 
through the lower layer 2. More precisely, one finds 
four rows of two inserts 7 for each, numbered 8 to 11, 
which extend respectively under the upper surface of 
the lower layer 2, under the second layer 3, under the 
layers 3 and 4 and under the three upper layers 3 to 5. 
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Generally, the inserts 7 comprise a portion in 
hydroxyapatite 12 and a portion in Plexiglas 13. The 
portions 12 and 13 have the same total height, but 
different respective heights in each row 8 to 11, such 
that the rays originating from the source 14 and which 
pass through each of the inserts 7, cross different 
combinations of thicknesses of the two materials. 
Furthermore, scattered radiation is produced analogous 
to that produced in a living being. This similitude is 
due to the constitution of the phantom 1 itself since 
the materials composing it, apart from simulating 
closely the soft tissues and the bones, have analogous 
proportions and distributions. In particular, the 
inserts 7 are sufficiently separated so that they do 
not receive scattered radiation coming from 
neighbouring inserts, but only from the base material 
of layers 2 to 5. 

Finally, this description will detail a method for 
estimation of scattered radiation, applicable for use 
in stages E2 and E5; however this method is part of 
another invention and is therefore only given in order 
to complete the description of this one and to prove 
the interest of the phantom 1. 

At the stage E2, radiation similar to that to be 
used on the subject and with flux 0 O is directed at the 
phantom 1, and a total radiation 0 t , the sum of a 
primary radiation 0 and the scattered radiation 0 d , is 
collected by the detectors placed behind the phantom 1, 
thus measuring the attenuation through it of the 
initial radiation 0 O . In particular, one obtains 
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measurements of radiation after passing across each of 
the inserts 7, both for high and low energy. 

The same measurements are repeated after setting 
absorbers in place, such as lead balls, between the 
source and the phantom in such a way as to intercept 
the primary radiation and only to measure the scattered 
radiation by the detectors behind the lead balls. By 
interpolation one then deduces the value of the 
diffused radiation 0 d in function of the heights of the 
two base materials crossed for each of the energies. 

The acquisition of the measurements through the 
subject also provides two ranges of values for total 
radiation 0 t at high and low energy. The radiation 0 d 
through the subject can be estimated according to the 
measurements behind the lead balls through the phantom 
1 and the relation between 0 t , 0d and 0 obtained through 
it for the two energies, assuming that the same 
relation applies through the subject. 

Although the invention is useful for conical 
beams, its use is not limited to these; one can 
conceive the stages of estimation and correction of 
scattered radiation being omitted in situations where 
it would be of less importance, since it would affect 
the measurements through the phantom 1 and the subject 
similarly, because of their similarity of structure, 
and its influence would disappear for the calculation 
of equivalent lengths Ai and A 2 . 



B 13769.3 JCI 



